Influence of energetically close orbitals on molecular high-order harmonic generation 
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We investigate the contributions from the 3(Tg and Itt^ and molecular orbitals in high-order 
harmonic generation in with particular emphasis on quantum-interference effects. We consider 
both the physical processes in which the electron is freed and returns to the same orbital, and those 
in which it is ionized off one orbital and recombines with the other. We show that the quantum- 
interference patterns observed in the high-order harmonic spectra are predominantly determined by 
the 3(jg orbital. This holds both for the situation in which only the Ittux orbital is considered, and 
the dynamics of the electron is restricted to the plane PxPz , or in the full three-dimensional case, if 
the azimuthal angle is integrated over and the degeneracy of Ittu is taken into account. 



I. INTRODUCTION 

In the past few years, high-order harmonic generation 
(HHG) has been extensively studied as a tool for at- 
tosccond imaging. In particular the possibility of bound- 
state reconstruction the attosecond probing of dy- 
namic processes in molecules and quantum interfer- 
ence effects Q has attracted a great deal of attention. 
This is a consequence of the fact that HHG is the result 
of the recombination of an electron, freed by tunneling 
or multiphoton ionization at an instant t' , with its par- 
ent molecule at a later instant t Since, in princi- 
ple, the electron may recombine with more than one cen- 
ter, high-harmonic emission at spatially separated sites 
takes place. Hence, information about the structure of 
the molecule in question is hidden in the HHG spectrum. 
In particular for diatomic molecules, this can be thought 
of as a microscopic counterpart of the double-slit exper- 
iment, in which maxima and minima arise due to the 
two-center interference Q. 

In many studies so far, it has been assumed that the 
electron is released from the hi ghes t-occupied molecular 
orbital (HOMO)[i, 0, i, i, 0, [O, [ll- This, however, 
has been disputed in recent investigations, in which it 
was shown that multielectron effects and the quantum 
interference of different ionization channels may play an 
important role [Tsj . Such effects may constitute a seri- 
ous obstacle towards ultra-fast molecular imaging. Apart 
from that, even if only the HOMO is considered as the 
initial state of the ionized electron, in many cases its de- 
generacy has a considerable influence on the HHG spectra 

In this paper, we investigate the influence of different 
molecular orbitals on the high-order harmonic spectra 
of diatomic nitrogen (iy2)- High-order harmonic gener- 
ation [^, 0, H, [HI, [ill and above-threshold ionization 
[3, [1, [3 i'^ -^2 have been extensively investigated in the 
literature, as, due to its large mass, its vibrational de- 
grees of freedom do not play a very important role and 
may be ignored to first approximation. In fact, it has 
been shown that, whereas for lighter species, vibration 
may lead to a considerable blurring of the two-center in- 
terference patterns, and a reduction in the high-harmonic 



or photoelectron yield, for molecular nitrogen such effects 
are not significant . 

Furthermore, in N2, the HOMO and the HOMO-1 or- 
bitals are energetically very close. This has several con- 
sequences. First, since they possess opposite parity, one 
expects a strong coupling between them. Second, since 
the tunneling probability is related to the bound-state 
energy, the processes in which the electron starts in the 
HOMO and in the HOMO-1 are comparable. Third, the 
electron may also leave from one orbital and recombine 
with the other, and, quantum mechanically, the transi- 
tion amplitude related to all physical processes involved 
will interfere. The influence of the HOMO-1 in the high- 
harmonic spectra of N2 has been recently observed 112;]. 

One should note, however, that, for N2, the HOMO 
and the HOMO-1 exhibit very distinct shapes and sym- 
metry. In fact, the former is a Scg orbital and the latter 
a ItTu orbital. Therefore, they are expected to behave 
differently as the alignment angle of the molecule with re- 
gard to the laser-field polarization is varied. Apart from 
that, the ItTk orbital is doubly degenerate. 

In our investigations, we employ the strong- field ap- 
proximation [20], and saddle-point methods. The tran- 
sition amplitudes obtained within this framework can be 
related to the classical orbits of an electron in a time- 
dependent field, and, yet, they retain information on 
the quantum interference between the possible physi- 
cal processes [2i|. Throughout, we employ the length 
gauge. Even though there is considerable debate about 
which gauge to employ, and the length gauge SFA leads 
to p otential-energy shifts whose meaning are not clear 
|la |. it has been recently shown that the two-center 
interference patterns are absent in SFA computations 
of the high-harmonic spectra using the velocity-gauge 

@,[il,[i9>]. 

This paper is organized as follows. In Sec. [Hi we 
provide the SFA transition amplitudes for the physical 
processes involved, for an exponential basis set involving 
Slater-type orbitals, and for a split- valence, gaussian ba- 
sis set. Subsequently, in Sec. IIIIl we compare the high- 
harmonic spectra obtained using both basis sets (Sec. 
nil Ap . and investigate quantum- interference effects be- 
tween the ScTg and l7r„ orbital (Sec. IIII Bp . Finally, in 
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Sec. IIVI we summarize the paper and state our main 
conclusions. 



II. TRANSITION AMPLITUDES 

Below we provide the HHG transition amplitudes, 
within the strong-field approximation. We base our ap- 
proach on the explicit expression in Ref . [20| , and employ 
atomic units throughout. 

The HHG amplitude is generalized to the case in which 
the active electron is initially in a coherent superposition 
of the 3CTg and the l7r„ orbitals. 

Explicitly, 

where the coefficients Cso-g, Ci7r„^ and Ci7r„^ give the 
weights of each state. One should note that the orbitals 
liTux and li^uy are degenerate, and possess the energy 
In the present model, we will neglect the processes 
in which the electron, immediately before ionization, is 
excited from the ItTu state to i<Jg, and, upon recombina- 
tion, decays from Scg to l7r„. 

Under these assumptions, the overall transition ampli- 
tude will be the sum 



M = ^ Afji, + c.c. 



(1) 



of nine terms. Explicitly, 
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x4i--)(p-f A(t'))exp[z533(t,i',p)], (10) 



where df"'\p) = (p| r.e, |3CTg) and di'""«^(p) = 
(p| r.e^ |l7r„^) , = x,y) are the components of the 
dipole matrix elements related to Bcr^ and l7r„ along the 
field-polarization axis. 

In the above-stated equations, one may distinguish two 
types of contributions. The amplitudes Mjj correspond 
to the processes in which the electron leaves a specific 
orbital, reaches a Volkov state |p-|-A(t')), propagates 
in the continuum from t' to t, and recombines from a 
Volkov state |p -I- A(t)) to the same orbital it left from. 
The amplitudes Mjv, j ^ v^ on the other hand, give the 
processes in which the electron leaves from one orbital 
and returns to another. The corresponding actions read 



S,,{t, i', p) = S(t, i', p) - E^(t - t'), 



(11) 



and 



S,,it,t',p) = S{t,t',p)-{E^t-E,3t'), j^v, (12) 



respectively, with 



sit,t',p)^nt-^j'dT [p + a(t)]' 



(13) 



My: 



dtjdt'j dv:^i-"^Hp+A(i)) 

—00 

x4i^-)(P + A(t')) exp[i522(t, t', p)], (6) 



In the above-stated equations, Ea, Ep refer to the bound- 
state energies, t to the recombination time, t' to the start 
time and p the intermediate momentum. For Sii,a = 
3ag, while for 6*22 and ^33, a = l7r„. For 5*12 (t, i', p), 
Ea — Esag and Ep = , whilst, for S2i{t,t',p), the 
situation is reversed, i.e., Ea — E^^^ and Ep = ^'So-g- 
Note that, due to the fact that the l7r„ orbitals are degen- 
erate, 5*12 = iSis, 5*21 = iS'31 and 6*32 = 5*22 = >5'33 = 6*23. 

We will compute the transition amplitude M employ- 
ing the stationary phase approximation, i.e., we will look 
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for values of t,t' and p that renders the actions in Eqs. 
(H])-® stationary. Apart from considerably simplify- 
ing the computations involved, this approach provides a 
physical interpretation of the amplitudes ilfjv in terms of 
electron trajectories. We compute the transition ampli- 
tudes employing a uniform saddle-point approximation. 
Details on the specific method used can be found in . 



A. Saddle-point equations 

Differentiating Su, S22 and S33 with respect to the 
ionization time t' and the recombination time t, we obtain 
the saddle-point equations 



and 



[p + Mt)Y 



n, 



(14) 



(15) 



where a — Sag for Sii{t,t' ,p) and a = l7r„ for S22 and 
533. Physically, Eq. pi)) gives the conservation of en- 
ergy at the instant of tunneling, and Eq. (jlSp expresses 
the fact that the electron recombines to the same state 
(either \3(Tg) or jlTTu)), releasing its kinetic energy upon 
return in form of a high-order harmonic of frequency fl. 
Finally, the condition dSjj/dp — yields 



dr [p + A(r)] = 



(16) 



Eq. ()16|) constrains the intermediate momentum of the 
electron, so that it returns to the site of its release. In 
the present model, this site is taken as the origin of our 
coordinate system, and is the geometric center of the 
molecule. Summarizing, the saddle-point equations (|14p - 
(fT6|) are related to the physical picture of an electron 
ionizing from either the HOMO or the HOMO-1 in N2 
and returning to the same state. 

The remaining actions Sji,{t,t' ,p), for j ^ v, lead to 
the saddle-point equations 



and 



[p+Mt)r 



Ea = ft, 



(17) 



(18) 



which indicate that the electron has left from one state 
and recombined with the other. For Si2{t,t' ,p) and 



S'i3(t,t',p), Ea = £'3^, and E, 
S2i{t,t',p) and^3i(i,t',P 
Ea = E^^^ and E/s = Es^^. 



■0 



^iTTuiWhile, for 



the situation is reversed, i.e.. 
For the remaining terms. 



TBI) remains the same 



Ea = Ep = E^^^^ so that Eqs. (|T4]) and are re- 
covered. Physically, this corresponds to the situation in 
which the electron leaves |l7r„2,) and returns to |lvr„j^) , or 
vice-versa. The return condition 
in this case. 



B. Orbital wavefunctions and dipole matrix 
elements 

Within the framework of the strong-field approxima- 
tion, all structural information about the molecule is em- 



bedded in the recombination prefactor d[ 
(p -f A(t)| r.e^ I*) , with ^ = 3ag or Itt^ 
space, this prefactor is given by 



p + A{t)) = 
In position 



1 



(27r)3/2 



d^rp • e^exp[-ip • r]'I'(r), (19) 



i.e., the component of idp'^{p) along the laser- field po- 
larization. In the following, we will construct the 
momentum-space wavefunction ^(p) for both orbitals. 
We will consider the linear combination of atomic or- 
bitals (LCAO) approximation and frozen nuclei. This 
implies that the position-space wavefunction reads 



«'(r) = ^7A4r + R/2) + (-l) 



i^air - R/2), 



(20) 

where R, la and rua denote the internuclear separation, 
the orbital and magnetic quantum numbers, respectively. 
For gerade and ungerade symmetry, Aa = \iTia\ and Xa = 
I ma I -f- 1, respectively. 

Throughout, we will use the length form of the dipole 
operator and neglect the terms growing linearly with the 
internuclear separation. Such terms are artifacts and 
come from the lack of orthogonality between the Volkov 
states and the field-free bound states. For a more com- 
plete discussion see [H, [13, [3 • 

The wave functions V'a(r) will be approximated by ei- 
ther exponentially decaying. Slater-type orbitals or by a 
gaussian basis set. In the former case. 



c.(2Ca)"°+^/^ ^„„ 
7(2^ 



-^»'^r""(0,0), (21) 



where Ua refers to the principal quantum number, and, 
in the latter, 



with 



V'fHr)=5]5a,^f)(r) 



ififir) = b^x^'^y'^^z^^ exp^C^^r^ 



(22) 



(23) 



The coefficients Cq, Caj and and the exponents <^ are 
extracted either from quantum chemistry codes, or from 
existing literature. 
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An exponential basis set has been recently employed 
in the literature 0, H, , while the use of gaussians is 
more widespread within the quantum chemistry commu- 
nity. In particular, a gaussian basis set exhibits several 
advantages. 

First, it allows an easier evaluation of the momentum- 
space wavefunction, which will be a central ingredient for 
computing the matrix elements di*-* (p + A(i)). Second, 
within the SFA framework, for exponentially decaying 
states, the ionization prefactor cfj^\p + A{t')) exhibits a 
singularity, according to the saddle-point equations (|14p 
and (|17p . In previous work, we have eliminated this sin- 
gularity by incorporating the prefactor d!"^\p + A{t')) 
in the action, and found out that it did not play a con- 
siderable role [2^. This singularity, however, is absent 
if gaussian wavefunctions are taken. Finally, in Hartree 
Fock computations there is an artifact that renders the 
ItTu orbital more loosely bound than 3(Jg. 

Explicitly, the Fourier transform of Eq. ([^D|) reads 



with 



r]{la,ma,p) = C+cos 



and 



p R 


-f iC- sin 


p R 






2 




2 


+ (-1; 







(25) 



(26) 



A generalized interference condition, which takes into 
account the structure of the orbitals in question, such 
as, for instance the s-p mixing in the Sag orbital, can 
be inferred from Eq. (|24p . Indeed, if we consider 
•ff — arctan(iC+/C_), then 



r]{la, ma, p) = ^JCl - Cl sin[^ + p • R/2]. (27) 

For ?y(/a: "TLa, p -l- A(t)), interference minima are present 
if 



I? + [p + A{t)] ■ R/2 = KTT, 



(28) 



where k denotes an integer number. This interference 
condition has been first derived in [1]. 

For Slater-type orbitals the individual wavefunctions 
are given by 



(~zp)'°2"°-'°Ca ^'"^'^'^ r(2 + la + ng) 
7(2^^ r(3/2-fna) 

X 2F1 (ai , q;2, as, "4)^^™° {Op, cj)p), (29) 



and the arguments of the hypergeometric functions read 
ai = \ + {la + na)l2, a2 ^ ai + 1/2, = la + 3/2, 04 = 
—p^/Ca- The angles are given by Op — cos~^ (p^/p) and 
(f)p = tan"i(py/p^). 



It is worth noticing that Eq. (|29|) is mainly employed 
in the description of a orbitals, since the spherical har- 
monics Y^°'{0p,4>p) are real for ma = 0. For tt orbitals, 
it makes physically more sense to employ real spherical 
harmonics, which are linear combinations of yj™" (0^, 0^) 
and Yf"^" {Op, 4>p). The explicit expressions for the real 
spherical harmonics are provided in Q. 

For a gaussian basis set the wavefunction V'a(p) reads 



with 



(30) 



(31) 



and k — x,y, z. Explicitly, 



*(P) =Z]^(^'^' + ^W)'^'^(P)' (24) Iwith 



MPk) = IC"^'' {x{Pk) + S(pfc)) , (32) 



x{Pk) = *(-i + (-i)^'=)p/cr(^a + ^) 



13- 



(33) 



and 



2(Pfc) = (1 + (-i)^--) v^r {(k) (^a; • (34) 

The arguments of the Hypergeometric functions are de- 
noted by 6 = (/3fe + l)/2 and ^ = -pl/{iCu)- In this 
work, we will be using p and s states, so that Eq. (I3ip 
will reduce to 

3 /2 

^.(P) = (-zpfe)'"^^575^exp[V/(4C.)]. (35) 
^ ^ (,1/ 

Therein, k = z, k — x and k — y for a, and tTj, 
states, respectively. The return condition (jl6p guarantees 
that the momentum p and the external field are coUinear. 
Hence, for a linearly polarized field Op is equal to the 
alignment angle Ol- 



III. HARMONIC SPECTRA 

In the following, we will present the high-order har- 
monic spectra. We choose the driving field as a linearly 
polarized monochromatic wave of frequency uj and am- 
plitude ojAq directed along the axis z. Hence, the corre- 
sponding vector potential is 



A{t) = Aq cos{Ljt)ez 



(36) 



For all situations, we consider starting times < i' < lot: 
confined to the first half cycle and the three shortest pairs 
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of orbits. For this particular field, using the saddle-point 
equation , the generalized interference condition 
may be expressed in terms of the harmonic order n as 



Eg 

LU 



LoR"^ cos^ 9l 



(37) 



where is the absolute value of the bound-state energy 
in question, n is an integer number, ^i, is the alignment 
angle, R is the internuclear distance and i) is defined in 
Eq. 



A. HOMO and HOMO-1 contributions 

In this section, we will make an assessment of the main 
differences encountered in the HHG spectra if the or- 
bitals are built employing a split valence, gaussian basis 
set, or exponentially decaying. Slater-type orbitals. For 
that purpose, we will concentrate on the recombination 
prefactor (ii*''(p -I- A(t)) and assume that the ionization 
prefactor (ii*''(p -I- A(t')) is constant and unitary. 

As a starting point, a direct comparison with the re- 
sults reported in 8] for the iag orbital will be performed. 
This orbital is known to exhibit a strong mixing between 
s and p states. Therefore, we will address the question of 
how such a mixing influences the overall interference pat- 
terns. In the context of the present article, this implies 
that we will consider the transition amplitude Mn, for 
which the electron leaves and recombines with the Zcjg 
orbital. 

In Fig.[Tl we display such results, either computed with 
a 6-3 IG gaussian basis set and coefficients obtained from 
GAMESS-UK m, or with Slater- type orbitals dH]), and 
the coefficients in [2^ [upper and lower panels, respec- 
tively]. The outcome of the split- valence computation. 



displayed in Fig. [l](a), exhibits a minimum which, for 
parallel molecular alignment, is near = 25w. This is a 
slightly higher harmonic order than that observed in [8] 
(see Fig. 4 therein). The minima observed for the indi- 
vidual s and p contributions, in contrast, agree with the 
results presented in Q (c.f. Fig.[TJ(b) and Fig.[TJ(c), re- 
spectively). This suggests that the s-p mixing possesses 
different weights in the present case and in Q . 

The spectra obtained with the Slater-type orbitals, on 
the other hand, are practically identical to the results 
in Q. This holds both for the minimum in the full 3(7g 
spectrum [Fig. [11(d)], which, for parallel alignment, is 
close to ri = 21a;, and for the patterns present in the 
s and p contributions [Fig. [T](e) and Fig. [TJ(f), respec- 
tively]. We have ruled out that this discrepancy is due 
to the slightly different ionization potentials employed in 
the two computations by performing a direct comparison 
for the same set of parameters (not shown). We have 
also found, employing GAMESS-UK and several types 
of basis sets, that the minimum at — 25a; is rather ro- 
bust with respect to small variations of £'30- and R. [2^ . 
Hence, in comparison to our computations, it seems that 




Alignment angle e[_ (degrees) Alignment angle 0[_ (degrees) Alignment angle e[_ (degrees) 



FIG. 1: High-order harmonic spectra for the HOMO in 
subject to a linearly polarized laser field of frequency 
LO — 0.057 a.u. and intensity / = 4x 10^*W/cm'^, as a function 
of the alignment angle 6l between the molecule and the field. 
The spectra in the upper panels have been constructed using 
a gaussian basis set and coefficients obtained from GAMESS- 
UK [2^ , while those in the lower panels have been built using 
Slater-type orbitals, and the coefficients in jS^]. From left to 
right, we display the spectra from the full Sctj orbital [panels 
(a) and (d)], the contributions from the s states [panels (b) 
and (e)], and those from the p states [panels (c) and (f)]. The 
bound-state energy of the HOMO and the equilibrium inter- 
nuclear distance have been taken from the respective com- 
putations. For the upper panels, -Eso-g = 0.63485797 a.u., 
while for the lower panels -Eao-g = 0.63495 a.u. In both cases, 
R = 2.068 a.u. 



the contributions of the s states to the spectra are slightly 
underestimated in 



24] 



In Fig. [21 we present the high-harmonic spectra com- 
puted assuming, instead, that the electron comes back 
and returns to the iTr^^: orbital, i.e., employing the tran- 
sition amplitude M22. The l7r„j, orbital should behave in 
a similar way and lead to the same spectrum, as it ex- 
hibits the same dependence with regard to the alignment 
angle. 

As in the 3crg case, we construct the bound-state wave- 
function either from Slater-type orbitals and the data in 
[23] or from a 6-31G basis set obtained from GAMESS- 
UK [2^. These results are displayed in Figs. [21(a) and 
[21(b), respectively. In both cases, we find that the two- 
center interference occurs at the very same harmonic or- 
der. Furthermore, apart from discrepancies in the overall 
intensity, the spectra exhibit a very similar substructure. 
Finally, in both cases, the yield drops considerably for 
parallel-aligned molecules. This is expected, as, if the 
angle Ol = 0, the it orbitals exhibit a nodal plane along 
the polarization axis. If the alignment angle increases, 
this nodal plane moves further and further away from 
the field-polarization axis, and the high-order harmonic 
yield increases. 
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Alignment angle e,^ (degrees) Alignment angle e^^ (degrees) 



FIG. 2: High-order harmonic spectra computed using the 
Ittux (HOMO-1) orbital in N2 as a function of the aUgnment 
angle Ol, for the same laser-field parameters as in the previ- 
ous figure. In panel (a), we considered a gaussian basis set 
and computed the coefficients with GAMESS-UK ^25], while 
in p anel (b) we took Slater-type orbitals and the data from 
[2J|. The bound-state energy of the HOMO and the equilib- 
rium internuclear distance have been taken from the respec- 
tive computations. For panel (a), -Ei7r„ = 0.61544 a.u., while 
for panel (b) Eit^^ = 0.65087981 a.u. In both cases, the in- 
ternuclear distance is 7? = 2.068. Note that in the Slater-type 
case the l7r„ orbitals are more loosely bound than the Sag 
orbital. 

B. Quantum interference of HOMO and HOMO-1 

We will now investigate which signatures the interfer- 
ence between the 3tTg and l7r„ leave on the high-order 
harmonic spectra. In all cases, we will consider both the 
recombination prefactor (ii*''(p -|- A(i)) and the ioniza- 
tion prefactor (ii*-'(p -|- A{t')). The ionization prefactor 
is important in this context due to the fact that an elec- 
tron reaching the continuum from a ct or a tt orbital be- 
haves in very different ways, with regard to the alignment 
angle 9l- In fact, for a orbitals, one expects ionization to 
be significant for small 9l and to be negligible for large 
values of this parameter. For tt orbitals, due to the pres- 
ence of the nodal plane, the opposite behavior is expected 
to occur. 



1. Two-dimensional model 

We will commence by addressing the situation for 
which (f)p = 0, i.e., we are restricting the dynamics of the 
problem to the PxPz plane. In this case, the initial wave- 
function is a superposition of the Sag and ItTux states 
only, i.e., Cit^^^ = 0. We consider that it is equally prob- 
able that the electron leaves from each of these states, 
i.e., C3a, = I/V2. 

In Fig. [31(a), we show the full spectrum, in which all 
the transition amplitudes Mjv, (j, v = 1,2) are summed 
coherently. Especially for small alignment angles, this 



spectrum exhibits a minimum very close to that obtained 
if only the iug state is taken. This minimum gets more 
and more blurred as the alignment angle 9^ increases. 
Possibly, this is the main influence of the 1t:ux orbital, as 
its contributions increase with 9^. In the following, we 
will investigate these patterns in more detail. For that 
purpose, we consider the quantum interference between 
specific processes. These results are depicted in the re- 
maining panels of Fig. [3l 

If the electron recombines with the ItTux orbital, re- 
gardless of where it started from [Fig. [31 (b)] , a very pro- 
nounced interference minimum is observed. This mini- 
mum occurs for the same harmonic orders as if only tt 
states are taken (c.f. Fig. [2]). This is expected, as high- 
order harmonic generation in the former case is only due 
to recombination with the li^ux orbital, even if two or- 
bitals are involved. 

Apart from that, the yield practically vanishes at 
9l — 0. This behavior is caused by the nodal plane 
which exists along the internuclear axis for the tTux or- 
bital. In this case, recombination for both the transition 
amplitudes M12 and M22, and ionization for the transi- 
tion amplitude M22, are strongly suppressed. For parallel 
alignment, this plane is along the laser-field polarization. 
As the alignment angle increases, this plane moves away 
from the field-polarization axis. Consequently, the yield 
increases. This explains why, in the overall spectrum, 
the minimum is determined by the iug orbital. For the 
parameters considered in this work, such a minimum lies 
mostly in the region of small alignment angles. 

In contrast, if only processes involving ionization from 
ItTux and recombination with 3(Tg, or vice- versa, are 
taken, the double-slit interference minimum is completely 
blurred [c.f. Fig. [31(c)]. This is due to the fact that 
both contributions are comparable, and exhibit minima 
for different harmonic orders. Furthermore, since either 
recombination with or ionization from a tt state is taking 
place, a strong suppression for 9l — is present. In Fig. 
[31 this is the only case for which we observed a complete 
disappearance of the double-slit minimum. Indeed, nei- 
ther for the processes involving only one state [31(d)], or 
starting at 3(Tg regardless of the end state [Fig. [31(e)] 
does the minimum completely vanish. However, a sharp 
minimum is only present if we take into account the pro- 
cesses in which the electron recombines with the same 
state. Concrete examples are Fig. [31(b), and Fig. [31(f), 
where the processes finishing at In^x and 3(Tg, respec- 
tively, are presented. In Fig. [31(f), we also notice that, 
for < 9l < 45°, the contributions from the a orbital 
are up to two orders of magnitude larger than those from 
the TT orbital [i.e.. Fig. [31(b)]. This is further evidence 
that the minimum is determined by the 3crg state. 



2. Three-dimensional case 

In a more realistic situation, one cannot restrict the 
electron dynamics only to the PxPz plane. In fact, there 
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FIG. 3: Contribution of different processes to the high- 
harmonic spectra, as functions of the alignment angle 9l, for 
the same field parameters in the previous figure. We chose 
</)p = so that the iTTuy orbital does not contribute. The 
dipole matrix elements have been computed using a gaus- 
sian basis set and GAMESS-UK [2^. In this case, -Eso-g ~ 
0.63485797 a.u, = 0.65087981 a.u. and R = 2.068 a.u. 

Panel (a): \Mii + M22 + M12 + M2if; panel (b): processes 
finishing at the Ittux orbital, i.e., IM21-I-M22P; panel (c): pro- 
cesses in which the electron starts at one orbital and recom- 
bines with the other, i.e., |Mi2 -I- M21P; panel (d): processes 
in which the electron starts from and returns to the same or- 
bital, i.e., |Mii -f M22P; panel (e): processes starting at the 
3(Tg orbital, i.e., \M2i + Muf; panel (f): processes finishing 
at the 30-9 orbital, i.e., \Mi2 + Mii|^. 



exist two TT orbitals which, even though they behave in 
the same way with respect to the alignment angle 9 l , are 
degenerate. Hence, they provide a completely different 
weight to the states ji/'o) from which the electron is re- 
leased and to which it returns. Furthermore, under many 
experimental conditions, the azimuthal angle (pp cannot 
be resolved. Thus, this parameter must be integrated 
over. 

Explicitly, the resulting spectrum is given by 



s{n) 



27r 



(38) 



TABLE I: Weights yV{a, /J,iy,j) for the contributions of the 
terms M^fjMj^ to the spectra, when integrated over the az- 
imuthal angle (j>p. 



above-stated sum consists of 81 terms. In general, the 
integrand in Eq ([55)) is of the form M*^Mjv. Its gen- 
eral dependence on the azimuthal angle (j)p is given by 
(sin (/>j,)''i (cos 0p)''^ , where the exponents ?7i,r?2 are inte- 
gers. Depending on such exponents, the contributions to 
the full harmonic spectrum carry different weights. For 
r]i or 772 odd, the contributions to the spectrum vanish. 
This implies that only the terms M*jM^^, for any j, i/, 
and M*jMj^, M*^Mj^, for j ^ survive. 

If the integral in ([55)) is carried out, one obtains 



For the specific problem addressed in this work, the 



^ W{a,P,u,])Mlp{p + A{t),6p) 

a,f3,u,j 

xMj,{p + A{t'),9p), (39) 

where M.jv{p + ^(t), 0p), t = t,t' is the transition am- 
plitude without the dependence on <pp. In Table 1, we 
provide it the weights yV(a, /3, j^, j) for each term in the 
sum psp. after integration over (j)p. 



In Fig. 131 we depict the high-order harmonic spec- 
tra obtained employing Eq. ([55)) . starting by the full 
spectrum [Fig. |3|(a)]. Therein, the minimum caused by 
the recombination of the electron with the Sag orbital is 
clearly visible, and the blurring due to the influence of 
the degenerate ItTu orbitals is even less pronounced than 
for its two-dimensional counterpart. At first sight, this 
is a counterintuitive finding, as, in the three-dimensional 
case, there are many more processes involving the latter 
orbitals. Possibly, this is a consequence of two main ef- 
fects. First, due to the presence of the nodal plane, for a 
broad range of alignment angles the contributions of the 
TT orbitals are strongly suppressed. Second, in general, 
the weights yV{a, (3,v,j) involving the Sag orbitals are 
larger than those involving the l7r„ orbitals only. This 
is, once more, a consequence of the geometry of the latter 
orbitals. 

In order to investigate this fact, we computed the 
high-order harmonic spectrum taking into account only 
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FIG. 4: Contribution of different processes to the high- 
harmonic spectra, as functions of the alignment angle 6l, 
for the same field parameters in the previous figure, tak- 
ing into account the degeneracy of the Ittu orbitals in a 
three-dimensional scenario, and non-resolved angle (jjp- The 
dipole matrix elements have been computed using a gaus- 
sian basis set and GAMESS-UK [2I]. In this case, £3^, = 
0.63485797 a.u, Ei„^ = 0.65087981 a.u. and R = 2.068 
a.u. Panel (a): all processes, i.e., the full sum in Eq. 1381 
panel (b): processes incorporating the Itvu orbitals only, i.e., 
\M22 + M23 + M32 + Mssf in l38l panel (c): processes in which 
the electron starts at any orbital and recombines with the Itvu 
orbitals, i.e.,[M22 -I- M23 -I- M32 -I- M33 -I- M21 + Msif; panel 
(d): processes in which the electron recombines with the Sag 
orbital, i.e., |Mii -I- M12 -I- Mi3|^ 



the latter contributions. Such results are displayed in 
Fig. 131(b). QuaUtatively, the spectrum obtained in this 
way is in perfect agreement with those displayed in Fig. [2l 
which have been computed using l7r„^ only, or with 
that shown in Fig. [31(b), which incorporates the pro- 
cesses in which the electron recombines at ItTux in a two- 
dimensional scenario. In fact, all such spectra exhibit a 
minimum above Q ~ 31w for 9p in the vicinity of zero, 
which moves towards the cutoff for Ol — 45°. The har- 
monic yield in Fig.|31(b), especially in the region of small 
alignment angles, is up to the three orders of magnitude 
weaker than the full spectrum. This huge discrepancy, 
however, would not cause any blurring in the full spec- 
trum. 

Potentially, the blurring may also be caused by the 
processes in which the electron is released from the Scg 
orbital and recombines with any of the l7r„ orbitals. In 
this latter case, a minimum near f2 = 31lu for 9^ in the 
vicinity of zero would also be present. Therefore, such 



processes must be incorporated. In Fig. [11(c), we con- 
sider the contributions from all possible processes finish- 
ing at the ItTu orbitals, regardless of where the electron 
left from. As expected, there is a substantial increase in 
the yield for small angles, in comparison to Fig. [41(b). 

Such increase is however not sufficient to match the 
contributions from the ScTg states to the full spectrum. 
In Fig. [51(d), we show that, for small alignment angles, 
the processes for which the electron recombines with the 
3(Tg state dominate. In fact, for < 9^ < 45° the yield 
in Fig. [41(d) is roughly one order of magnitude larger 
than that displayed in Fig. [11(c). This is not obvious, as 
there are twice as many more processes contributing to 
the yield in this latter case, namely six against three. A 
direct comparison of Figs. [41(a) and [41(d) also shows the 
above-mentioned dominance for small angles. For larger 
angles, the contributions from the l7r„ orbitals start to 
play a more significant role and there is an increase in the 
blurring. In all cases for which the electron only starts 
from or recombines with the l7r„ orbitals [Figs. [11(b) and 
[11(c)], there is a strong suppression of the yield for par- 
allel alignment. This is due to the fact that the nodal 
plane along the molecular axis coincides with the laser- 
polarization axis in this case. 



IV. CONCLUSIONS 

We considered the influence of two closely lying molec- 
ular orbitals on the high-order harmonic spectrum from 
N2: the ItTu and Sag orbitals. We employed a very sim- 
ple model, in which the strong-field approximation has 
been modified in order to incorporate the situations in 
which an electron leaves from the Sug orbital and recom- 
bines with ItTu and vice-versa. We have also included 
the degeneracy of the l7r„ orbital. We made a detailed 
assessment of the contributions of all possible processes 
to the high-order harmonic spectra. 

The main conclusion to be drawn from this work is 
that the shape and the two-center interference patterns 
observed for the high-order harmonic spectra from 7V2 
are mainly determined by the Sag orbital, even though 
the ItTu orbitals are energetically very close. The main 
effect of the latter orbitals is to introduce some blurring 
in the interference minimum determined by Sag. 

Physically, this is due to the particular geometry of the 
ItTu orbitals. Indeed, for small alignment angles 9l, these 
orbitals exhibit a nodal plane close to the polarization 
axis, so that tunneling and recombination are strongly 
suppressed. Hence, in this region, the high-order har- 
monic spectra are mainly dominated by the Sag orbital. 
We have verified that this dominance extends up to ap- 
proximately 9l = 45°. For the parameters considered in 
this paper, the two-center minimum occurs within this 
region, so that it is mainly determined by the Sag state. 

Furthermore, due to their nontrivial dependence on the 
azimuthal angle, the l7r„ orbitals carry less weight when 
this parameter is integrated over. Interestingly, even if 
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a three-dimensional computation is carried out and the 
degeneracy of the l7r„ orbitals is considered, this angular 
dependence outweighs the fact that there are more pro- 
cesses in which the electron recombines with one of the 
l7r„ orbitals. 

We have also shown that, due to the above-mentioned 
non-trivial angular dependence, the influence of such or- 
bitals is over-estimated if the dynamics of the problem 
is reduced to the PxPz plane, i.e., if the azimuthal angle 
(j}p is chosen to be vanishing. Such an approximation has 
been extensively used in the literature (see, e.g., [1| in 
with HHG from the tt^ orbital of the O2 molecule has 
been computed). This is not an obvious result, as two- 
dimensional models do not consider the degeneracy of the 
Ittu orbital. 

Finally, it is worth mentioning that the findings of this 
paper agree qualitatively with recent results obtained em- 
ploying more sophisticated methods, such as Dyson or- 
bitals and many-body perturbation theory . Therein, 
it has been shown that many-electron effects did not play 
a significant role in the bound-state reconstruction of N2 , 
and that the information retrieved from the spectra was 
mostly related to the Sug orbital. 



It may, however, be possible to identify the influence of 
the ItTu orbitals by looking at effects for perpendicular- 
aligned molecules, or relatively large alignment angles. 
In this case, the contributions from Sag to the high-order 
harmonic spectrum are not expected to obfuscate those 
from l7r„. In fact, recently, the influence of the latter 
orbitals on the HHG spectrum of N2 has been identi- 
fied experimentally for perpendicular-aligned molecules, 
in form of a maximum at the rotational half-revival . 
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